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Abstract 

The fluctuation determinant which determines the preexponential factor of the 
transition rate for minimal bubbles is computed for the electroweak theory with 
sin Qyy = 0. As the basic action we use the three-dimensional high-temperature 
action including, besides temperature dependent masses, the T$ 3 one-loop con- 
tribution whcih makes the phase transition first order. The results show that this 
term (which has then to be subtracted from the exact result) gives the dominant 
contribution to the one-loop effective action. The remaining correction is of the 
order of, but in general larger than the critical bubble action. The results for 
the Higgs field fluctuations are compared with those of an approximate heat ker- 
nel computation of Kripfganz et al., good agreement is found for small bubbles, 
strong deviations for large thin-wall bubbles. 
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1 INTRODUCTION 

The electroweak phase transition is at present the object of extensive investiga- 
tions 0. If the phase transition is first order, which is possibly the case if the 
mass of the Higgs boson is not too large, the phase transition occurs via bub- 
ble nucleation. Bubble nucleation can have various consequences for cosmology 
in the early universe. The possibility of baryogenesis in bubble walls has been 
investigated recently by many authors (see e.g. |2|, H), reheating after the phase 
transition could be mediated by bubble nucleation and subsequent coalescence, 
the creation of inhomogeneities by bubble formation could be observable (see e. 
g. ||, |5|, U for representative discussions of the physics of bubble nucleation and 
growth) . 

Bubble nucleation is described usually within the reaction rate theory formu- 
lation of Langer or, equivalently, the semiclassical approach to quantum field 
theory by Coleman and Callan || ||. This formulation requires the existence 
of a saddle point in configuration space, the minimal bubble, with one unstable 
mode, possible zero modes and real frequency fluctuation modes. The leading 
term in the tunnelling rate is given by the negative exponential of the minimal 
bubble action, the corrections arise from integrating out the fluctuations in the 
Gaussian approximation, leading to a fluctuation determinant prefactor whose 
negative logarithm is the 1-loop effective action. If the leading approximation 
is good this prefactor should be of order 1, substantial prefactors have however 



been found in the case of the sphaleron transition, both from bosonic [10, 11, 12 



as also fermionic [13|] fluctuations. It is therefore of interest to investigate how 



strongly these prefactors modify the leading order approximation to the bubble 
nucleation rate. 

Here we present an exact computation of the bosonic fluctuation determinant 
of the critical bubble. As the basic action is determined by the usual Higgs poten- 
tial with just one minimum at the classical expectation value, some fluctuation 
effects have to be included already at the tree level in order to allow for mini- 
mal bubble solutions. The exact fluctuation determinant should then reproduce 
those in order to justify this modification of the leading order action. Following 
the basic work of Coleman and Weinberg such modified actions have been 



proposed by many authors || [L5], [H| [T7| and used to describe the bubble nucle- 
ation in leading order. To be specific we use here the one given by Dine et al.|| 
which was also the basis of a recent approximate computation of the 1-loop Higgs 
fluctuations by Kripfganz et al. |18| . 



The plan of this paper is as follows: In the next section we will introduce the 
model and set up the basic relations for the bubble nucleation rate. In section 3 
we will discuss the structure of the fluctuation operator, in particular its partial 
wave decomposition. The computation of its determinant, based on a very useful 
theorem, will be described in section 4. In the final section we will present some 
results and conclusions. 
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2 Basic relations 



The three-dimensional high-temperature action is given, in the formulation by 
Dine et al. M, by 
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Here the coordinates and fields have been rescaled as 
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The vacuum expectation value v(T) is defined as 

- 2 (T) = ^(T 2 -T 2 ). 

T is the temperature at which the high-temperature potential Vht changes its 
extremum at $ = from a minimum at T > T to a maximum at T < T . The 
temperature dependent coupling of the three-dimensional theory is defined as 

gT 
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In terms of the zero temperature parameters we have mw = gvo/2, mu = V 2At>o 
with vq = 246 GeV and we use the definitions of Dine et al. || modified by 
setting Qw = and therefore mw = rnz '■ 
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In terms of these parameters the high-temperature potential is given by 
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The rescaling Eq. ( |2.2| ) with the scale v(T) makes sense only for T < T . On 
the other hand the high-temperature potential has, before rescaling, a secondary 
minimum at = v(T) with 
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This minimum is degenerate with the one at $ = at a temperature defined 
implicitly by 



T c = To/yjl - Ey\ Tc . (2.9) 
Tc marks the onset of bubble formation by thermal barrier transition. In the 



work of Hellmund et al. and Kripfganz et al. [|18j the vacuum expectation 
value of the broken symmetry phase v(T) is chosen for the rescaling of the fields, 
i.e. in Eq. ( |2.2|) v (T) is replaced by v(T), the high-temperature coupling constant 
Eq.(|2.4j) is redefined analoguously and denoted^] as g$(T). By this change of scale 
the high-temperature potential changes as well; it becomes [] J2T 



V ht (&&) = ^ (V<fr) 2 - e(T)($t$)3/2 + ( 3 £(T) _ 2) ^^ (2 1Q) 



with 



The action and its rescaling differ slightly from that of Hellmund et al. PU 



and of Kripfganz et al. ||18|| . In contrast to the former we do not mimic the 



influence of a Debye mass by decoupling the longitudinal degrees of freedom. In 
contrast to the second one we include only the $ 3 contribution of gauge field 
and would-be Goldstone degrees of freedom as in Ref. [16]. This form of the <3> 3 



contribution was found to yield a good approximation for the exact results in the 
case of the sphaleron |10], |llj, at least in the case m H /mw <C 1. We will find, 
indeed, that this term dominates the effective action. 

The process of bubble nucleation is - within the approach of Langer [[jj and 
Coleman and Callan || |{J, followed by the work of Affleck |2l|], Linde and 
others - described by the rate 

V I V = ^ f eM~S)J- 1/2 . (2.12) 

Z7T \ Z71 J 



Here S is the high-temperature action, Eq. (2.1), with the new rescaling, min- 
imized by a classical minimal bubble configuration (see below), J is the fluctu- 
ation determinant which describes the next-to-leading part of the semiclassical 
approach and which will be defined below; its logarithm is related to the 1-loop 
effective action by 



1 
2 

Finally is the absolute value of the unstable mode frequency. 



Sljj = ~ In J. (2.13) 



2 Our notation differs from the one of Refs. jl8], g0| 
3 We do not introduce a tilde for the rescaled fields. 
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The classical bubble configuration is described by a vanishing gauge field and 
a real spherically symmetric Higgs field $(r) = |3>|(V) which is a solution of the 
Euler-Lagrange equation 

- d>"(r) - -$'(r) + = (2.14) 

with the boundary conditions 

lim $(r) = and $'(0) = . (2.15) 

r— »oo 

This differential equation can be solved numerically e.g. by the shooting method. 
The solution will be denoted as H (r). 



3 Fluctuation analysis 

In terms of the action S the fluctuation operator is defined generally as 

8 2 S 

■^ ab = KA \<t>=<t>bubble' (3.16) 

d<fi a d(p b 

where 4> a stands for the various gauge and Higgs field components and 4>bubbie is 
the field configuration of the minimal bubble. An analoguous derivative, taken at 
= (pvac = defines the vacuum fluctuation operator M°. In both configurations 
the gauge fields vanish, the Higgs field is given by 



* = H (r) ; (3.17) 




in the bubble configuration and vanishes in the vacuum. 

The fluctuation determinant J appearing in the rate formula is defined by 

J ~^M«- (3 ' 18) 

Here the symbol det" denotes the determinant with removed translation zero 
modes and with the unstable mode frequency replaced by its absolute value. 

The analysis of fluctuations of the minimal bubble can be related to a similar 
analysis performed recently for the electroweak sphaleron without gauge fixing 



in Ref. |23[ and in the 't Hooft-Feynman background gauge in Ref. |TT|. We 
will use this latter analysis. One can take over the fluctuation operator with two 
modifications which represent at the same time essential simplifications: 

- the high-temperature effective potential has to be modified from the one in Eq. 
( P77D to the one in Eq. ( gTTJ) ; 

- the sphaleron and the broken symmetry vacuum configurations are replaced by 
the bubble and the symmetric vacuum configurations defined above. 
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Furthermore we use here (see Eq. ( |2.2| )) for the coordinates the scale (gv) 1 

2/gv used in Ref. 



instead of the scale M w x 
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The expansion of gauge and Higgs fields reads then [ 1 ! 
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(3.19) 



Here the fields denoted with small letters , a^h and 4> a are the fluctuating fields. 

Before we discuss fluctuations we have to fix the gauge. We work here in the 
't Hooft-Feynman background gauge. The gauge conditions read 



1 
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The total gauge-fixed action S t is obtained from the high-temperature action by 
adding to it the gauge-fixing action 



Sgf ~ gi(T) J d3x 2 :Fa:Fa 



(3.21) 



and the Fadeev-Popov action 
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It is the action S t = Sht + S g f + Spp which has to be used in the definition of the 
fluctuation operator ( |3.16| ). 

The Hilbert space of fluctuations decomposes into subspaces defined by the 
symmetries of the background field. The fluctuation operators given below have 
been derived from those of Ref. ||11|| . This analysis was based on a K spin basis 
(K = J + I). Alternatively one might have used here simply an analysis based on 
ordinary spin, i.e. an expansion where the Higgs field, the Fadeev Popov field and 
the time components of the gauge fields are expanded with respect to spherical 
harmonics and the space components of the gauge fields with respect to vector 
spherical harmonics xY; m ,rV^ m and LYj" 1 . 

The electric components of the gauge field and the isovector (would-be Gold- 
stone) components of the Higgs field form a coupled (3 x 3) system. The fluctu- 
ation operator can be written in the form M = M° + V. The free operator M° 
is diagonal. It consists of free partial-wave Klein-Gordon operators 



d 2 



2 d , l n (l n + 1) 



dr 2 r dr 



+ 



(3.23) 



with masses m n given by (0, 0, m#) respectively for the three components and 
with centrifugal barriers corresponding to angular momenta l n given analoguously 
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by (/ + 1, 1 — 1, 1). The nonvanishing components of the potential are 

V u = V 22 = H 2 /A 
V 33 = H 2 /4+( \ T /4g 2 ) (4H 2 -3eH ) 
I I + 1 dH 

= v - = -1WTTlf (3 - 24) 



v ~ = WaTTlT- (3 ' 25) 

For Z = the second component is absent due to the vanishing of the vector 
spherical harmonic tVYq. These amplitudes have a triple degeneracy due to 
isospin besides the ordinary degeneracy (21 + 1) from spin. 

The fluctuation operator for the scalar part of the Higgs field is given by 



d 2 2d 1(1 + 1) 



M = -— - -— + " v " - '' + m\ + V M (r 



dr 2 r dr r 2 H 
V M = ^(l2H 2 -6eH ) (3.26) 



™H = i^( 3e " 4 )- 

This channel being an isosinglet its degeneracy is just (21 + 1). 

The time components of the gauge fields, the Fadeev Popov fields and the 
magnetic components of the vector potentials all satisfy the same equation 

M^s = w 2 ^5 • (3-27) 



It consists of a free massless partial wave Klein Gordon operator 

d 2 _ 2 d 1(1 + 1) 

dr 2 r dr r 

and a potential 



V 55 = 3 (3.29) 

which vanishes exponentially as r — > oo. There is no I = component of the 
magnetic vector potential since the vector spherical harmonic LY™ vanishes. In 
the fluctuation determinant all of these contributions cancel, only the s-wave 
Fadeev- Popov contribution survives, due to the lack of its magnetic counterpart. 
It is triply degenerate due to isospin and has to be subtracted. 

The partial-wave decomposition of the fluctuation operator decomposes also 
its determinant, 

J = J2( 2l + 1 )Ji- (3-30) 
We now need a method for computing numerically the determinants of the partial 
wave fluctuation operators. Such a method has been developed recently by V. G. 
Kiselev and the author and will be presented briefly in the following section. 
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4 The fluctuation determinant of the electroweak 
bubble 

A very fast method for computing fluctuation determinants is based on a theorem 
on functional determinants; references to earlier work and an elegant proof are 
given in Ref. p4]. Generalized to a coupled {n x n) system it can be stated in 



the following way: 

Let f(z/, r) and f°(V, r) denote the (n x n) matrices formed by n linearly 
independent solutions f^(v,r) and f"°(v,r) of 

(M tJ + v 2 )f«(v,r)=0 (4.1) 

and 

(M°. + ^ 2 )/j>,r) = 0, (4.2) 

respectively, with regular boundary conditions at r = 0. The lower index denotes 
the n components, the different solutions are labelled by the greek upper index. 
Let these solutions be normalized such that 

limf(i/,r)(f°(z/,r))- 1 = 1 . (4.3) 

r— >0 

Then the following equality holds: 

= det(M + ^) = detffor) 
1 } det(M° + i/ 2 ) ™detf°(z/,r) 1 ' ; 

where the determinants on the left hand side are determinants in functional space, 
those on the right hand side are ordinary determinants of the n x n matrices 
defined above. If the theorem is applied at v — it yields the desired ratio 
of fluctuation determinants J = J7(0). The consideration of finite values of 
v is necessary in the discussion of zero modes. The theorem has been applied 
for computing the 1-loop effective action of a single scalar field on a bubble 
background in Ref. and of a fermion system on a similar background in Ref. 
2lJ previously. It was found to yield very precise results, in addition to providing 



a very fast computational method. 

In the numerical application the solutions /° were written as [27 



/;w-K+^:(r)kM (4.5) 

with the boundary condition h"(r) — > as r — > 0. The values l n and K n = 



m\ + z/ 2 depend on the channel as specified in the previous section. This way 
one generates a set of linearly independent solutions which near r = behave 
like the free solution as required by the theorem which then takes the form 

3{v)= hmlndet{^ + ^(r)}. (4.6) 
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The functions h"(r) satisfy the differential equation |57 



^-Mr) +(- + ^nP^4) U<* n {r) = V nn/ (r) + K,{r)) (4-7) 



dr* nw \r n ii n (K n r) dr nX ' ™ v v * ■ • ( x 



which can also easily be used for generating the functions order by order in 
V. In particular, if this differential equation is truncated by leaving out the term 
h®> on the right hand side, one generates the first order contribution to h which 



is the tadpole term. For more technical details we refer to Refs. [Pq , |2q . 

With the partial wave fluctuation operators given in the previous section the 
application of the theorem to the case of the electroweak bubble is straightfor- 
ward. Some points to be considered are 

- the subtraction of the divergent tadpole graphs 

- double counting of gauge and would-be Goldstone fluctuations 

- removing the translation zero mode 

- removing a particular gauge zero mode. 

We will discuss these briefly. We will add also some remarks on details of the 
numerical computation. 

4.1 Tadpole diagrams 

The high-temperature three-dimensional theory has only linear divergences of 
the form of tadpole diagrams which renormalize the mass term of the Higgs 
field. They have to be subtracted in the numerical computation to obtain finite 
results. This was done in each partial wave, for which the tadpole contribution 
may be computed p5| , [26[j either by solving a truncated differential equation or 



27J as an analytic expression using the partial wave Green function. After these 
contributions have been subtracted, the partial wave contributions converge as 
l// 2 and have a finite sum. 

Of course this contribution has to be added back, after having been regularized 
and renormalized. Part of these diagrams have already been taken into account 
in the renormalization of the four dimensional theory and in giving the vacuum 
expectation value ( |2.3| ) of the Higgs field a quadratic temperature dependence. 
Some terms linear in the temperature survive however and contribute |10[ [II], 
[28f| (after dividing by the temperature) to the 1-loop effective action, i.e. the 
logarithm of the fluctuation determinant. If the mass of the field in the loop 
is rrii and its coupling to the external field is described by the potential Vi their 
contribution to the effective action is given by — rrii/ 871 J d 3 xVi(r). The fluctuating 
gauge fields have vanishing mass and do not contribute. However we receive 
contributions from the fluctuating Higgs fields. The mass circulating in the loop 
is then uih which is - including the temperature dependence and rescaling - given 
by equation ( |3.26| ). The potentials are V33 with triple isospin degeneracy and V44. 



S 



So we have to restitute the terms 



S)Z l tad = -^fJ drr- (jif 2 + 3^(4tf 2 - 3eH ) ] (4. 
for the gauge fields and 



StL = / drr%^- 2 {2Hl - eH ) (4.9) 



for the Higgs field to the 1-loop effective action. 

We should like to remark on two slight inconsistencies of this procedure. The 
first one concerns our choice of the high temperature action. We have adopted the 
action of Dine et al. M since it sets a certain standard and since it has been used 
also in Ref . |18| to which we want to compare part of our results (and which shares 
the inconsistency). In this action the T 2 term does not include the contribution 
of the Higgs loop tadpole. This can be seen from the coefficient D in Eqs. ( |2.5| ) 
which should include a term m 2 H in addition to Zm^ v + 2m 2 . The contribution was 
neglected already in Ref. |17| "taking the Higgs boson sufficiently light". Since 



the expression is dominated by the top quark contribution whose fluctuations are 
not included at all here this omission may be tolerated at the present level of 
accuracy. In a more refined analysis it should and can easily be remedied. 

The second point is the fact that we have included already 1-loop effects into 
the tree level action, so that part of our computation is now at the 2-loop level, 
without constituting a complete and systematic 2-loop analysis. This applies 
in particular also to the tadpole terms for which this could be a more severe 
problem since they are the finite remnants of divergent graphs. We can appeal 
here only to an argument - common in many perturbative calculations - that 
possible inconsistencies are of higher order and acceptable at an intermediate 
level as they will be cured in a complete higher order analysis. 



4.2 Double counting of gauge field fluctuations 

As mentioned in section 2 we are working with an action that contains already the 
part of the 1-loop effective potential induced by integrating out the gauge field 
and would-be Goldstone boson fluctuations. These are present in the temperature 
scale factors and couplings and appear especially in the high temperature effective 
potential as the term proportional to $ 3 . While the T 2 contribution to the 
vacuum expectation value ( |2.3| ) comes from the tadpole diagrams and has been 
taken into account along with these, the $ 3 term is contained in our exact 1- 
loop effective action. In order to avoid double counting it, this term has to be 
subtracted from our numerical results. The incorporation of this term into the 
tree level action was necessary in order to obtain a first order phase transition 
and bubble solution. If this was a good leading order approximation the gauge 
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field action should be well approximated by this term. This is indeed the case (see 
below) but this implies also that the remaining gauge and would-be Goldstone 
field contributions are small differences of large terms, and that they cannot 
therefore be expected to be very precise. 



4.3 Translation zero mode and unstable mode 

Translation invariance is broken by the classical solution, so a zero mode appears. 
It occurs in the I = 1 partial wave of the fluctuation operator of the isoscalar 



part of the Higgs field. It is easily removed using the prescription given in |p5 
one applies the theorem mentioned above at finite v and defines 

Ji=x,m 99 s = lim lim In (^^\) (4.10) 



where k = \jv 2 + m 2 H . Removing three eigenvalues u 2 = gives the fluctua- 
tion determinant J the dimension (energy)^ 6 . The rate gets then a dimension 
(energy) 3 = 1 / (length) 3 . An additional dimension energy = 1/time comes from 
the unstable mode prefactor (see ( |2.12|) ).The numerical computation is based on 
energy units gv(T) which are given in the Tables below. 

The unstable mode makes the determinant of the p-wave contribution neg- 
ative. Replacing it by its absolute value means just to revert the sign of the 
determinant before taking the logarithm. We note that, in contrast to Ref. 



and in analogy to Refs. |T0| pj| , we do not remove the zero mode from the 
fluctuation determinant. 



4.4 Gauge zero mode 

Though we have imposed a gauge condition there is one residual gauge degree of 
freedom. It is analoguous to a constant gauge function for the free theory. Indeed 
in the latter case a constant gauge potential A(x) = go does not contribute to the 
vector potential and is therefore not eliminated by the gauge condition d^a^ = 0. 
In the case of the bubble background field there is a similar but nontrivial mode 
which satisfies the background gauge condition and is therefore not eliminated 
by it. It manifests itself as a zero mode in the electric system for 1 = 0. The 
form of this mode (and the fact that it is really an exact zero mode) was found 
after extended numerical experiments. It is given by a gauge function g(r) which 
satisfies the same differential equation as the electric and Fadeev Popov modes, 
Eq. CT , i-e. 

g"+ 2 ~g'-^9 = o. (4-11) 
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With regular boundary condition at r = g(r) becomes constant as r — > oo, in 
analogy to the free case. Then the functions 

^i(r) = -2g\r) 

Vs(r) = H g(r) (4.12) 

satisfy the coupled system for the electric modes at I = which is given explicitly 
by 



2 „/. _ AT ( ATJ2 



^' + -^-m^ 3 = ^3 + ^(4# 2 - 3eifo)^3 - • (4-13) 
r 4 4g z 



It can be checked easily that this gauge zero mode satisfies the background gauge 
condition ( |3.20|) and is therefore not eliminated by it. Since this zero mode is not 



due to a symmetry broken by the classical solution as the translation mode it 
cannot be handled in the usual way. On the other hand we observe that precisely 
for the s-wave the Fadeev-Popov contribution has survived; furthermore to each 
Fadeev-Popov mode with finite energy, i. e. a solution of 

€ + -^-^5 = -^5 (4.14) 
r 4 

there is a solution of the electric s-wave system constructed exactly as that for the 
gauge zero mode, i.e. Eq. ( |4.12|) with g replaced by ip 5 . So there is a cancellation 
of all electric modes of this type with the corresponding Fadeev-Popov ones, 
except for the mode with uj^ = 0. There is of course a solution of the Fadeev- 
Popov equation at this energy, but it is "singular" at infinity, going to a constant 
there. The corresponding mode in the electric system is normalizable, however, 
since only its derivative is involved in ipi and its product with the exponentially 
decreasing function Hq in if) 3. The cancellation between the s-wave electric modes 
(|4.12[) and the Fadeev-Popov ones can be extended therefore to the zero mode if 
the boundary condition at r — > 00 for the latter ones is replaced by ip' 5 (r) — > 0. 
This can be done in analogy with the procedure described in the previous section 
by computing the fluctuation determinant of the Fadeev-Popov mode at finite v 
via 

JwpM = In lim ( ^T^rl • (4-15) 

Then the Fadeev-Popov system at I = exhibits a zero mode as well, the limit 

lim(j; = o,e^) - Ju*M v )) (4-16) 

is finite and defines the s-wave part of the fluctuation determinant. In this way 
the Fadeev-Popov term cancels all unwanted longitudinal electric modes for I = 
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- including the one with frequency zero. We note that the change of boundary 
condition as r — > oo affects only the s-wave and only for massless fields. The 
definition ( |4.15|) yields results identical to the usual one ( f4.4| ) if I ^ and/or the 
fields are massive. 



4.5 Some numerical details 

The analysis was performed as described in previous publications |25|, |26|. Con- 
tributions of angular momenta up to l max = 30 were computed numerically, the 
higher ones were included by performing a power fit Al~ 2 + Bl~ 3 + C7 -4 through 
the last 5 computed contributions and by adding a corresponding sum from l max 
to oo. This was done already at lower values of Z, treating the highest included 
angular momentum as the actual value of l m ax- The resulting expressions were 
found to be independent of I within typically four significant digits for / > 20. 

A more subtle point is the extrapolation to r = oo implied in Eq. (}4.4|) . In the 
previous analyses [EH EB1 the fields had finite mass and the approach to r = oo 



was exponential. For the massless fields the Bessel functions ii(fvr) are replaced by 
r / (21 + 1)!! and the functions and approach their asymptotic value only as 
hoo + const ./ r . The extrapolation was performed using this Ansatz. An exception 
occurs in the electric p-wave system, where /if picks up a logarithmic dependence 
on r due to the cross term with h% (r) on the r.h.s. of Eq. ( |4.13j ) which decreases 
only as 1/r 2 . However this logarithmic dependence being strictly proportional 
to #3 + /^(oo), i.e. to the third row of the matrix, it does not contribute in the 
determinant, as also observed numerically. 

The tadpole contributions were computed in two ways, once by solving a trun- 
cated differential equation as described in p5[ |26|] and performing the analoguous 
extrapolation, and once as an integral using the partial wave Green function. In 
comparing the two results the extrapolation was found - for the tadpole contri- 
butions - to be reliable to four significant digits typically. 

Judging the accuracy of the results from the stability with respect to varying 
extrapolations as r — > oo and for large I we would think that the purely numerical 
part is accurate to 1 %. The restituted tadpole contributions are given by the 
expressions (|4.8| ) and ( f4.9|) whose evaluation implies simple numerical integrals, 
they can be considered as exact analytic expressions. This restitution implies 
no delicate cancellations. However even with a precision of 1% for the numeri- 
cal results the final values of the gauge field contribution have substantial errors 
since the numerical part plus the tadpole contribution is almost cancelled by the 
analytic $ 3 contribution. Unfortunately, in contrast to the sphaleron computa- 
tion [PH| , the cancellation ist not merely one between two analytic expressions 
- the tadpole and $ 3 contributions - but between the numerical results and the 
analytic $ 3 contribution. 



12 



5 Results and conclusions 



The numerical results are given in Tables 1 to 4. Here Table 1 is based on the 
values Higgs and gauge boson masses m# = mw = 80.2 GeV, and a value of 
the top mass of m t = 170 GeV. For the vacuum expectation of the Higgs field 
we used v = 246 GeV and for the gauge coupling the value g = .6516. For 
the computation of Table 2 the values mg = 60 GeV, m t = 170 GeV were used. 
Table 3 corresponds to values mj? = 60 GeV and m t = 140 GeV, Table 4 to values 
rriH = 80.2 GeV and mt = 140 Gev; these latter Tables are presented in order 
to compare with results obtained in Ref. |18| using the heat kernel expansion. 
The values for the temperature chosen correspond to 10 equidistant steps of the 
quantity e(T), defined in Eq. ( |2.11| ), between the onset of bubble nucleation at 
e = 2 and the critical temperature Tq where bubble nucleation ends at e = 4/3. 



This choice is equivalent to the choice of Kripfganz et al. |T8[ who parametrize 
this range of temperatures by a variable y taking values between and 1. Since 
Kripfganz et al. use a somewhat different effective potential, the relation between 
y and e is not precise, it is essentially given by y = 3 — 2e which we use as a 
definition of 'our' y. At small y the bubbles are large with thin walls, for y — 1 
the bubbles are small and have thick walls. 

Tables 1 and 2 are split into a part 'a' which contains the essential parameters 
for the minimal bubble and, in the last column, the nucleation rate R without 
fluctuation corrections. The part 'b' contains the fluctuation corrections, i.e. the 
1-loop effective action. The results for m t = 170 GeV are given separately for the 
isoscalar part of the Higgs field as Sl~ l and for the system of would-be Goldstone 
fields and gauge fields ('gauge field contribution' for short) as Sg~ l , respectively. 
We give also separately the parts which were obtained by the numerical analysis 
described in section 4, denoted as S] l j num and SgZ l num , respectively. The difference 
between and S]^} num is the tadpole contibution S]^} tad of Eq. (|4.9|), and 
analoguously for the gauge field. Note that the tadpole contributions to the 
Higgs field are substantial. The gauge field contribution Sg~~ l contains the $ 3 
part discussed in the previous section. The numerical value of this term is given 
in the column '$ 3 '. This term should be close to the gauge field contribution, and 
it is indeed. So the basic action used for computing the bubble profiles represents 
a reasonable approximation to the exact 1-loop effective action. The gauge field 
contribution has to be reduced by this term since it would be double-counted 
otherwise. The net gauge field contribution is denoted as AS 1 ^ and given in 
the last column. The correction to the rate can be simply obtained as a factor 
exp(— ASljj) where ASljj = + AS^~'. The dimension energy 3 , here in 
units of gv, is already included in the minimal bubble rate R. One sees that the 
fluctuations lead to a substantial suppression of the nucleation rate. 

While the result that the effective action for the gauge fields is well approx- 
imated by the effective potential, i. e. the $ 3 term, is very rewarding a less 
comfortable feature appears if one compares the 1-loop effective action with the 
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tree level action S. If the saddle point approximation which forms the basis of 
transition rate formula fl2.12|) is justified, then the 1-loop action should be smaller 
than the tree level one. This is not the case. Large 1-loop corrections were found 
already by Kripfganz et al. Jl8| when computing the 1-loop effective action for 
the Higgs-field fluctuation only. We compare our results to theirs in Tables 3 
and 4. Since these authors define the fluctuation determinant differently - they 
remove the unstable mode - we give, besides our result S}f iggs , the expression 
In(v4/T 4 ) where A is the square root J~ x l 2 of the fluctuation determinant with 
translation and unstable modes removed. The results are close to each other for 
small e or y ~ 1, i .e. for small thick- wall bubbles. For small y, i.e. for large 
thin-wall bubbles, our exact results are systematically larger than the approx- 



imate ones of Ref. ||18|| . The question of finding reliable analytic estimates is 
certainly an interesting one, especially the order in which the terms of the heat 
kernel expansion are summed. In |26[] it was found that a summation by the 



number of derivatives ("derivative expansion") yields very precise results if the 
mass of the fluctuation is much larger than the inverse size of the background 



field configuration. In Ref. |18| the terms are summed with respect to powers of 



the heat kernel time. The deviation at small y could be due to the fact that large 
bubbles with thin walls have a very substantial derivative contribution. It will be 
interesting to make a more systematic analysis of various analytical approaches. 
The comparison of the Higgs effective action with the leading minimal bubble 
action is less favorable than found in Ref. [P8]. This is even more the case if the 



gauge loops are included, as one sees from the previous Tables 1 and 2. 
In conclusion we state three essential features of our results: 

- The 1-loop effective action is substantial, of the order of and larger than the 
leading order minimal bubble action. This sheds some doubt on the applicability 
of the semiclassical transition rate theory. 

- The sign of the 1-loop effective action is such that the transition rate is sup- 
pressed. 

- The 1 -loop "$ 3 " contribution which has been incorporated into the basic effec- 
tive potential is reproduced rather well by the 1-loop action. This means that this 
term in the effective potential describes relevant features of the effective action. 

It will be interesting to pursue this subject further; it could be of interest to 
try a selfconsistent extremalization of the sum of leading and 1-loop action. This 
would certainly reduce the total suppression. Furthermore it will be interesting 
to see how the inclusion of the fermion determinant affects the transition rate. 
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Table Captions 



Table la Parameters of the minimal bubbles for mu = 60 GeV and m t = 170 
GeV. The results are given as a function of temperature in equidistant steps of 
the variable e (Eq. ( |2.11| )). v(T) is the temperature dependent vacuum expecta- 
tion value of Eq. ( |2.3|) , Aj- the temperature dependent renormalized $ 4 coupling. 
S is the minimal bubble action (or energy divided by T). uo 2 _ is the square of 
the frequency of the unstable mode, given in units of gv(T) 2 . The last column 
contains the logarithm of the nucleation rate without the 1-loop corrections. 

Table lb S]^} num is the 1-loop of the isoscalar part of the Higs field as obtained 
in the numerical analysis. Sl~ l is the total Higgs part of the 1-loop effective ac- 
tion, obtained from S^-num by adding the tadpole contribution Slz l ta d- SgZ l nurn is 
the 1-loop gauge and would-be Goldstone field action obtained by the numerical 
analysis, Sg~ l is again obtained by including the tadpole contribution. The next 
colum gives the $ 3 term as included into the high temperature action. AS^ - ' 
is the gauge field action after subtraction of this $ 3 contribution. ASljj is the 
total effective action after removing the $ 3 contribution. 

Table 2a The same as Table la for mu = rn w = 80.2 GeV and m t = 170 
GeV. 

Table 2b The same as Table lb for mu = my/ = 80.2 GeV and m t = 170 
GeV. 

Table 3 Comparison of the Higgs field effective action with approximate re- 
sults by Kripfganz et al. for mu = 60 GeV amd m t = 140 GeV. The first entries 
are as defined in the previous Tables. The quantity A is the Higgs part of the 
fluctuation determinant with removed unstable mode. Our results are compared 
to the one of Ref. [|I8"| , marked with the subscript KLS. 

Table 4 The same as Table 3 for m H = 80.2 GeV and m t = 140 GeV. 
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Tables 



T [GeV] 


e 


y 




A T x 10 2 


S 


uj 2 _ x 10 3 


ln^GeV 4 ]) 


94.557 


1.933 


.1 


48.82 


3.309 


1114.2 


-0.1947 


-1098.7 


94.529 


1.866 


.2 


50.53 


3.310 


278.47 


-0.7911 


-264.21 


94.495 


1.800 


.3 


52.37 


3.311 


121.40 


-1.819 


-107.83 


94.455 


1.733 


.4 


54.33 


3.312 


65.218 


-3.338 


-52.13 


94.405 


1.663 


.5 


56.47 


3.314 


37.875 


-5.438 


-25.20 


94.347 


1.600 


.6 


58.72 


3.316 


22.686 


-7.814 


-10.45 


94.276 


1.533 


.7 


61.19 


3.319 


12.891 


-9.760 


-1.223 


94.191 


1.466 


.8 


63.85 


3.322 


6.3659 


-9.928 


4.423 


94.089 


1.400 


.9 


66.75 


3.325 


2.1002 


-7.007 


7.028 



Table la 



y 


ql-l 


Ql-l 


ql-l 
g—num 


ql-l 


$3 


AS 1 ' 1 

^ gauge 




.1 


-389.2 


505.0 


-85974 


-102494 


-104089 


1595 


2100 


.2 


-81.21 


133.6 


-8457 


-9952 


-10170 


218 


351.6 


.3 


-28.77 


64.25 


-1984 


-2285.5 


-2338.3 


52.83 


117.1 


.4 


-11.78 


28.30 


-663.6 


-739.33 


-757.91 


18.59 


46.89 


.5 


-4.41 


28.85 


-263.44 


-277.53 


-287.56 


10.03 


38.88 


.6 


-7.54 


23.04 


-118.53 


-113.26 


-121.51 


8.25 


31.29 


.7 


1.48 


19.64 


-54.93 


-42.63 


-51.13 


8.50 


28.14 


.8 


3.11 


17.71 


-24.85 


-10.04 


-19.51 


9.47 


27.18 


.9 


4.79 


16.96 


-10.04 


5.53 


-5.14 


10.67 


27.63 



Table lb 
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T [GeV] 


e 


y 


dlGeV] 


A T x 10 2 


S 


uj 2 _ x 10 3 


ln(i?[GeV 4 ]) 

\ L J / 


115.725 


1.933 


.1 


39.75 


4.973 


601.83 


-0.294 


-587.87 


115.702 


1.866 


.2 


41.16 


4.974 


151.04 


-1.190 


-138.32 


115.675 


1.800 


.3 


42.67 


4.975 


65.944 


-2.733 


-53.91 


115.642 


1.733 


.4 


44.29 


4.976 


35.314 


-5.030 


-23.76 


115.602 


1.663 


.5 


46.05 


4.977 


20.560 


-8.173 


- 9.42 


115.555 


1.600 


.6 


47.91 


4.978 


12.329 


-11.73 


- 1.61 


115.498 


1.533 


.7 


49.96 


4.980 


7.066 


-14.65 


3.09 


115.428 


1.466 


.8 


52.19 


4.982 


3.450 


-14.88 


5.82 


115.345 


1.400 


.9 


54.49 


4.984 


1.145 


-10.51 


6.64 



Table 2a 



y 


ql-l 
h—num 


Ql-l 


ql-l 
g—num 




$3 


AS* 1 "' 

gauge 




.1 


-389.27 


501.83 


-44642 


-54007 


-56066 


2060 


2561 


.2 


-81.76 


132.93 


-4463.4 


-5236.8 


-5513.5 


276.7 


409.6 


.3 


-29.40 


63.66 


-1067.5 


-1194.4 


-1270.6 


76.21 


139.9 


.4 


-12.33 


39.39 


-362.64 


-377.30 


-409.36 


32.06 


71.45 


.5 


-5.01 


28.22 


-148.57 


-137.31 


-155.97 


18.66 


46.88 


.6 


-1.36 


22.42 


-69.30 


-52.00 


-66.02 


14.04 


36.46 


.7 


0.87 


19.03 


-33.64 


-15.43 


-27.78 


12.35 


31.38 


.8 


2.51 


17.09 


-16.35 


1.26 


-10.56 


11.82 


28.91 


.9 


4.20 


16.36 


-7.70 


8.91 


- 2.80 


11.71 


28.07 



Table 2b 
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T [GeV] 


e 


S 


u 2 _ x 10 3 


v[GeV] 


ql-l 
Higgs 




V 1 / KLb 


97.311 


1.933 


1336.1 


-0.1686 


56.88 


505.2 


-513.40 


-115. 


97.268 


1.866 


336.57 


-0.6965 


58.87 


134.2 


-141.56 


-82.0 


97.216 


1.800 


146.77 


-1.601 


61.00 


64.60 


-71.40 


-50.6 


97.153 


1.733 


78.407 


-2.955 


63.30 


40.20 


-46.54 


-36.7 


97.078 


1.663 


45.935 


-4.777 


65.75 


29.13 


-35.06 


-29.5 


96.988 


1.600 


27.934 


-6.891 


68.38 


23.23 


-28.83 


-25.4 


96.879 


1.533 


15.554 


-8.613 


71.24 


19.83 


-25.15 


-23.0 


96.748 


1.466 


7.6707 


-8.758 


74.34 


17.89 


-23.01 


-21.8 


96.590 


1.400 


2.5103 


-6.155 


77.70 


17.18 


-22.30 


-22.0 



Table 3 



T 


e 


S 


x 10 3 


v [GeV] 


ql-l 
° Higgs 


ln(£) 


In (&) 

\ T J KLS 


123.815 


1.933 


598.02 


-0.296 


43.12 


490.0 


-500.0 


-89. 


123.783 


1.866 


153.96 


-1.175 


44.63 


132.83 


-142.0 


-81.4 


123.744 


1.800 


67.333 


-2.695 


46.26 


63.712 


-72.32 


-52.0 


123.697 


1.733 


36.144 


-4.949 


48.01 


39.485 


-47.61 


-38.3 


123.640 


1.666 


21.072 


-8.035 


49.92 


28.306 


-36.06 


-31.0 


123.572 


1.600 


12.563 


-11.58 


51.95 


22.428 


-29.84 


-26.8 


123.490 


1.533 


7.1533 


-14.44 


54.16 


19.047 


-26.18 


-24.4 


123.391 


1.466 


3.5336 


-14.67 


56.57 


17.118 


-24.06 


-23.1 


123.270 


1.400 


1.1142 


-10.25 


59.22 


16.363 


-23.30 


-23.3 



Table 4 
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